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The Gamma-ray Burst Monitor (GBM) on board the Fermi Gamma-ray Spaee 
Teleseope has triggered on over 300 Terrestrial Gamma-ray Flashes (TGFs) sinee 
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its launch in June 2008. With 14 detectors, GBM collects on average ~ 100 counts 
per triggered TGF, enabling unprecedented studies of the time profiles of TGFs. 
Here we present the first rigorous analysis of the temporal properties of a large 
sample of TGFs (278), including the distributions of the rise and fall time of 
the individual pulses and their durations. A variety of time profiles are observed 
with 19% of TGFs having multiple pulses separated in time and 31 clear cases 
of partially overlapping pulses. The effect of instrumental dead time and pulse 
pile-up on the temporal properties are also presented. As the observed gamma- 
ray pulse structure is representative of the electron flux at the source, TGF pulse 
parameters are critical to distinguish between relativistic feedback discharge and 
lightning leader models. We show that at least 67% of TGFs at satellite altitudes 
are significantly asymmetric. For the asymmetric pulses, the rise times are al¬ 
most always shorter than the fall times. Those which are not are consistent with 
statistical fluctuations. The median rise time for asymmetric pulses is ~3 times 
shorter than for symmetric pulses while their fall times are comparable. The asym¬ 
metric shapes observed are consistent with the Relativistic Feedback Discharge 
(RFD) model when Compton scattering of photons between the source and Fermi 
is included and instrumental effects are taken into account. 
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1. Introduction 

Terrestrial Gamma-ray Flashes (TGFs) were first observed by the Burst and Transient Source 
Experiment (BATSE) on board the Compton Gamma-Ray Observatory [Fishman et al, 1994] 
as brief, intense flashes of y-ray emission coming from the Earth’s atmosphere and are linked to 
electrical activity in atmospheric regions in or above thunderstorms (e.g., Fishman et al. [1994]; 
Smith et al. [2008]). The leading model for TGE production relies on the relativistic runaway 
electron avalanche (RREA) process [Gurevich et al., 1992; Roussel-Dupre and Gurevich, 1996; 
Lehtinen et al, 1999], in which energetic electrons are accelerated in the presence of the strong 
electric field of a thunderstorm and produce an avalanche of relativistic electrons, which emit 
y-rays via bremsstrahlung radiation. Positrons and photons created by the runaway avalanche 
may in turn contribute to producing more energetic electrons in a process known as relativistic 
feedback, which can more easily explain the observed fluxes [Dwyer, 2003, 2007]. Eor a recent 
review of TGEs see Dwyer et al. [2012]. 

Since their discovery, the properties of these events have been studied with BATSE [Fishman 
et al., 1994], the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) [e.g., 
Grefenstette et al., 2009], the Astro-rivelatore Gamma a Immagini EEggero (AGIEE) [e.g., 
Marisaldi et al., 2010] and the Gamma-ray Burst Monitor (GBM) and Earge Area Telescope 
(EAT) on board the Fermi Gamma-ray Space Telescope [e.g., Briggs et al., 2010; Grove et al, 
2012]. These observations have shown TGEs to be characterized by very short durations and 
hard spectra. They have generally been found to consist of one or in some cases multiple pulses 
of y-ray emission, with each pulse having a typical duration of less than 1 ms. BATSE detected 
76 TGEs over 9 years which exhibited varied temporal behaviour, including several TGEs with 
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multiple distinct pulses. The TGF durations ranged from sub-millisecond to several millisec¬ 
onds [Nemiroff et al, 1997; Gjesteland et al, 2010]. However the triggering timescale of 64 ms 
used for BATSE caused an observational bias towards multi-pulse and longer events [Grefen- 
stette et al, 2008; Gjesteland et al, 2010]. The observed BATSE temporal properties were also 
shown to have been heavily modified by instrumental effects [Grefenstette et al, 2008; Gjeste¬ 
land et al, 2010]. RHESSI amassed a large sample of 820 TGEs between 2002 and 2008 and 
shows a much larger proportion of short single-pulse events than BATSE, with almost all hav¬ 
ing durations < 1 ms [Grefenstette et al, 2009]. A new analysis of AGIEE observations shows 
a median duration of ~370//s [Marisaldi et al, 2013]. The comparison of the results from dif¬ 
ferent missions is complicated by the intrinsic differences in the detectors used and the various 
methods used to determine the duration. 

Despite these TGE observations, a rigorous analysis of the temporal properties of individ¬ 
ual TGE pulses has not yet been performed. This is due to the low counting statistics and/or 
low temporal resolutions of the above instruments which do not allow such an analysis to be 
performed. Instrumental effects such as dead time and pulse pile-up also affect the observed 
temporal properties, further complicating such studies. 

The time profile of the gamma ray emission at the source follows the runaway electron 
flux. Dwyer [2012] suggests that relativistic feedback naturally makes roughly symmetric TGE 
pulses, regardless of the pulse duration, and asymmetric pulses may be observed due to Comp¬ 
ton scattering as the photons propagate to spacecraft altitude. It has been suggested that TGEs 
may be produced by mechanisms associated with the stepping of lightning leaders [Moss et al, 
2006; Celestin and Pasko, 2011]. In such a case, short TGE durations (~50 /us) may be consistent 
with the temporal dispersion associated with the Compton scattering of photons produced by an 
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instantaneous TGF source, and longer TGFs may be due to overlapping emissions. Constrain¬ 
ing the shapes and durations of TGF pulses is therefore of utmost importance in distinguishing 
between TGF production models. 

GBM has triggered on hundreds of TGFs since its launch in June 2008. Despite being op¬ 
timized for y-ray burst observations, GBM is particularly well-suited to TGF studies. The 
large bismuth germanate (BGO) detectors are particularly important for the detection of the 
higher-energy photons associated with TGFs and enable the accumulation of a large number of 
counts. On average, the 14 GBM detectors register ~100 counts per triggered TGF. Combined 
with an absolute timing accuracy of several yus, this allows the temporal properties of TGFs to 
be constrained with unprecedented accuracy. Briggs et al. [2010] found that the pulses in a 
small sample of 12 GBM TGFs (including at least one unrecognised Terrestrial Electron Beam 
(TEB)) tend to have symmetric or asymmetric shapes and are well-fit with Gaussian and log¬ 
normal functions. Fishman et al. [2011] presented a qualitative overview of the time profiles of 
the first 50 TGEs triggered by GBM but did not perform a quantitative analysis of individual 
pulse shapes. 

In this paper we perform a rigorous quantitative temporal analysis of TGEs by fitting functions 
to a large sample of individual TGE pulses. We present the rise times, fall times and durations 
derived from these pulse fits and quantify pulse asymmetry. We also investigate the impact of 
instrumental effects on the observed temporal properties of TGE pulses. The sample described 
in this paper comprises 278 of the first 300 TGEs which triggered the GBM instrument, covering 
the time period from the launch of Fermi in June 2008 up to 06 November 2012. In addition 
to those TGEs detected by GBM via in-orbit triggering, an offline search for TGEs which did 
not trigger the instrument has been implemented [Briggs et al, 2013]. The TGEs in this new 
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sample are typically much less fluent than the triggered sample (see Figure 4d in Briggs et al. 
[2013]). The counting statistics required to constrain the fit parameters are not present in most 
of the untriggered events and we therefore limit our analysis to the sample of triggered GRBs. 

Section 2 describes the GBM instrument. In Section 3, the pulse-fitting procedure is de¬ 
scribed. The results are presented in Section 4 and discussed in Section 5. 

2. Gamma-ray Burst Monitor 

The Fermi Gamma-ray Space Telescope was launched in June 2008 and has two instruments 
on board, the Large Area Telescope (LAT) dedicated to high-energy y-ray observations from 
20 MeV to > 300 GeV [Atwood et al, 2009] and the Gamma-ray Burst Monitor (GBM). GBM is 
comprised of 12 sodium iodide (Nal) scintillation detectors and two BGO scintillation detectors, 
covering the energy ranges ~ 8 keV to 1 MeV and ~ 200 keV to 40 MeV respectively [Meegan 
et al, 2009]. The Fermi satellite has a nearly circular orbit of 565 km altitude and an inclination 
of 25.6 ° which means that for a large fraction of its orbit it is located above the tropics, where 
thunderstorms are prevalent. 

GBM has an on-board triggering process, based on a count rate increase which is statistically 
significant in two or more detectors within a given timescale and over a particular energy range. 
In the event of a trigger, high-resolution time-tagged event (TTE) data are telemetered to the 
ground covering the time period from ~ 30 s before the trigger time to ~ 300 s after the trigger 
time. These data consist of individual photon events in each detector, tagged with their arrival 
times at 2 /rs resolution (with an absolute accuracy of several //s) and associated energy. The 
rate of TGF triggers is currently ~ 1 every 3-4 days. Due to their short durations and extremely 
high fluxes, the high photon rates of TGFs detected by GBM result in significant dead time 
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losses {Briggs et al, 2010]. The degree to which dead time affects our results is discussed in 
Section 3.1. 


3. Pulse-fitting Analysis 

Previous studies of a small sample indicate that GBM TGF pulses tend to be either symmetric 
or have faster rise times than fall times and are well described by simple Gaussian or lognormal 
functions [Briggs et al, 2010]. To investigate the true pulse properties of a large TGF sample, 
a function with greater flexibility is fit to the individual pulses of 278 GBM TGFs, in addition 
to a constant background rate. The Norris function [Norris et al, 1996] is chosen as it allows 
for a wide variety of pulse shapes, including those with longer rise times than fall times. It is 
formulated as: 


m 


^max exp 
^max £Xp 




, t <tp 

, t>t„ 


( 1 ) 


where Rmax is the maximum rate in counts per second per detector, tp is the peak time (s), cr,. is 
the rise time constant for t < tp (s), CTf is the fall time constant for t > tp (s) and v is a measure 
of the sharpness of the pulse, where lower values indicate a more peaked pulse. A peakedness 
value of 2 represents a two-sided Gaussian function (the cr parameter in the Norris function is 
related to the cr of a Gaussian by a factor of two in the exponent). As it was found that the full 
set of parameters was typically poorly constrained for weaker TGFs in the sample, v is frozen 
to 2 for all Norris fits. The best-fit parameters are those which maximise the Poisson likelihood 
function. A complete description of the analysis process is given in Briggs et al [2010]. 

The rise (tr) and fall (t/) times of the pulses are derived from the fit parameters. For each 
pulse, tr and tf are defined as the time intervals between the 10% and 90% flux levels of the 
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rising and falling edges of that pulse respectively, and are given by Eqns. 2 and 3: 


(Tr 


LV 








log 


1 

oT 


V 






yj 


log 


0.9 




(2) 

( 3 ) 


LV ' ' / V ' ' n 

For V = 2, the rise/fall times are related to the respective cr by a multiplicative factor of 1.21. 
Errors on these pulse properties are obtained by propagating the uncertainties on the function 
parameters. To ensure robust statistics, it is necessary to perform the pulse-fitting analysis over 
the full energy band of the Nal and BGO detectors. As there is known spectral evolution in 
TGFs, with a general trend for events to soften over time which has been attributed to Compton 
Scattering (e.g., Grefenstette et al. [2008]; 0stgaard et al. [2008]), the pulse parameters will 
vary depending on the energy interval studied. The derived rise time should be generally speak¬ 
ing, independent of the energy band used, as the rise time interval is dominated by unscattered 
photons. For the fall time, the energy band is more important, as a lower energy interval should 
have a longer fall due to the greater proportion of scattered events. Therefore, we stress that the 
rise times reported here should be representative of the incident rise time, whereas the fall time 
will be slightly longer than the photon fall time above 1 MeV and slightly shorter than the fall 
time for < 1 MeV. 

The TGFs described in this paper are those which triggered the GBM instrument and all 
events have full TTE data coverage. The pulse fits are performed using the combined TTE data 
of all 12 Nal and 2 BGO detectors over the full GBM energy range (8keV to 40 MeV) at the 
intrinsic 2/rs resolution of the data without dead time or pulse pile-up correction. The Fermi 
GBM data are public and are available at the Fermi Science Support Center^. A web catalog of 
GBM TGF lightcurves is also available^. 
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The pulses were also fit with normal (symmetric) and lognormal (asymmetric) functions as 
defined in Briggs et al. [2010]. While the more flexible Norris function allows for a single 
function to be fit to all TGF pulses, the normal and lognormal fits are also presented here as 
they have a more physical basis that can be more easily interpreted than the Norris function. 
It also allows the present results to be compared with previous studies and the results using 
simpler functions may be useful as the basis for further work or TGF simulations. The criterion 
for preferring a normal or lognormal fit to the data is based on the Akaike Information Criterion 
(AIC, Akaike [1974]), which allows for the comparison of the goodness of fit of statistical 
models. It is defined as: 


A/C = 2k-2 log X, (4) 

where k is the number of free parameters in the model and X is the maximized value of the 
likelihood function for that model. The model which minimizes the AIC is considered to be 
the best fit to the data. Therefore if the maximum likelihood is identical for two models, the 
AIC selects the one defined with the smaller number of parameters. The AIC is also used to 
determine if an additional pulse is required by the fit for those TGFs with possible multiple 
overlapping pulses. 

3.1. Dead time and Pulse Pile-up 

As the data we fit are subject to dead time and pulse pile-up losses {Chaplin et al, 2012], the 
peaks of the of the TGF pulses presented here are underestimated. Using the BGO detectors, 
Tierney et al. [2013] performed corrections for dead time, pulse pile-up and detection efficiency 
in a model independent manner for 106 TGFs detected by GBM. They concluded that the pri¬ 
mary observational effect was on the observed maximum rate and fluence, and that temporal 
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properties were not significantly affected. However, that work did not use the Nal detectors, 
for which the effect of pile-up is greater due to the higher proportion of overflow counts in the 
Nal data compared to the BGO data (overflow counts, which are above the maximum digitised 
energy, incur a dead time of 10/i 5). 

In order to quantify the degree to which the observed parameters of TGFs detected by GBM 
are influenced by instrumental effects, we have performed Monte Carlo simulations. Eight rep¬ 
resentative TGF pulse shapes (4 symmetric and 4 asymmetric) of varying duration and intensity 
were each simulated 1000 times. The resulting synthetic TGFs were then passed through a filter 
which emulates the detection process and includes the effects of dead time and pulse pile-up. 
The pulses were then fit with the Norris function and the observed pulse parameters compared 
to the input. A full description of the simulation process and analysis is given in Appendix A. 

As expected, the observed peak rate was heavily modified for both symmetric and asymmet¬ 
ric TGFs, with the observed value ranging from ~ 10 - 50 % of the input depending on both the 
duration and initial rate. For the temporal parameters, the mean values of the measured distri¬ 
butions are broadly consistent with the input values, with a tendency to shift to longer times. 
This shift is generally small, with the input value usually contained within +1 cr of the mean 
of the dead time filtered distribution. The exception is the fall times for the asymmetric TGFs, 
which show the same trend but to a greater degree. This elongation has the effect of artificially 
increasing the observed duration. However, the effect is small, even for very asymmetric high- 
rate TGFs, with the mean of the observed duration distribution typically differing from the input 
value by 20-30 ixs. These results are in agreement with those of Tierney et al. [2013]. 

4. Results 
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The T 90 dunitioii mecisiire hus previously been used for TGF pulses [e.^., ct al,^ 2010]}, 
whieh is the time it takes to observe 90% of the eounts in a pulse, starting and ending when 5% 
and 95% of the total eounts have been observed. This measure is highly dependent on an 
aeeurate estimation of the baekground level and is impossible when there is pulse eonfusion for 
overlapping pulses. Additionally, as the Tgo is modified by dead time, the observed Tgo should 
be eonsidered an upper limit on the true duration. Based on the Norris pulse fits presented in 
this paper, a new TGF pulse duration measure is proposed, defined as the 5% to 95% level 
integration of the pulse fit. Figure 1 shows the 7^90 value as a function of the duration derived 
from the pulse fit. The proposed duration measure agrees with the Tgo values measured for 
isolated pulses for the 226 TGF pulses for which both duration measures were available. An 
advantage of this new measure is that the effect of dead time can be estimated from simulations. 
The simulations presented in Appendix A show that for the highest-rate TGF simulated, the 
effect of dead time is to increase the duration by ~ 10 %. 

Of the first 300 TGF-like events detected in orbit by the GBM flight software, 22 are single- 
pulsed events with durations derived from the Norris pulse fits of longer than 1 ms. Of these, 8 
are confirmed TEBs with an additional very likely TEB. Since the underlying physical processes 
producing TEBs are thought to be understood [Dwyer et al, 2008; Briggs et al, 2011], we 
exclude these 9 events from the subsequent pulse-fitting analysis. The weaker long events are 
more difficult to classify; many are likely to be TEBs. In order to avoid contaminating the 
sample with unidentified TEBs, these 13 events are also excluded from the analysis. 

Of the remaining 278 TGEs, 55 have multiple emission episodes, consisting of two or more 
pulses. The maximum number of statistically significant pulses observed in a single TGE is 
three, with nine TGEs having three resolvable pulses. There are 31 instances of overlapping 
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pulses which are statistically significant. Of course, we cannot exclude the possibility of other 
multi-pulse TGFs in the sample for which additional pulses cannot be resolved. An example of 
a TGF with multiple and overlapping pulses fit with the Norris function is shown in Figure 2. 

The pulse-fitting analysis results in 348 individual pulses, to which two cuts were applied. 
For overlapping pulses the fit parameters are correlated and ill-constrained. A cut removing all 
of these was applied, leaving 280 isolated pulses. The second removed pulses which had poorly 
constrained rise and fall times (defined as uncertainties greater than 50% of the value). This 
resulted in a total of 250 isolated pulses with well constrained parameters. 

A scatter plot of rise time against fall time obtained from the Norris fits is plotted in Figure 3 
for all of the pulses that passed the cuts. In this and all subsequent figures, the error bars 
indicate the 68% confidence intervals. In general, the rise time is shorter than or comparable to 
the fall time for TGF pulses. A similar plot is shown in Figure 4 for pulses fit with Gaussian 
and lognormal functions. Consequently the sample is divided into those which are consistent 
with symmetry (33%; red on Figures 3 and 5), and those which are asymmetric (67%). Here 
we define an event as consistent with symmetry if the interval corresponding to the rise time + 
its uncertainty overlaps with that defined by the fall time + its uncertainty. Otherwise, the event 
is considered to be asymmetric. 

Sixty-three percent of the total sample have shorter rise times than fall times (blue), with 
just 4% of the sample having longer rise than fall times (green). The distributions of the rise 
times and fall times derived from the parameters of the pulse fits are shown in Figure 5. It is 
evident that the rise time distributions for the symmetric and asymmetric pulses with tr < tf are 
considerably different, with median values of 116/iS and 43 yus, respectively. The fall times are 
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comparable for symmetric and asymmetric pulses, having median values of 122/is and 120yus, 
respectively. 

Just ten pulses have longer rise times than fall times (green on Figures 3 and 5). Most of these 
are consistent with statistical fluctuations from the symmetric group as discussed in Appendix 
A3. One or two may be due to unrecognized overlapping pulses that were fit with a single Norris 
pulse. 

The temporal variability of TGF pulses derived from the pulse parameters can also be used 
to constrain the source emission radius [Dwyer, 2008]. A source with characteristic dimension 
R cannot vary faster than t ~ R/c due to the time required for a ‘turn-on’ signal to propagate 
and turn on 100% of the source [Carroll and Ostlie, 2006]. Various models predict different 
propagation velocities at the source and therefore different source radii. Using the rise time as 
the shortest temporal variation of a pulse, an upper limit can be placed on the source radius. 
Table 1 gives the derived upper limits on the source radii for the three pulses with the shortest 
rise times for the speed of light and for the typical propagation velocity of a lightning stepped 
leader. For a lightning stepped leader, there are two velocities to consider, the speed at which the 
lightning extends while it is advancing and the average velocity of the extension of the lightning 
which includes the many pauses and steps. In this analysis, we are interested in the former, 
which we take as ~ 1 x 10^ m/s [Rakov and Uman, 2003, p. 135]. The fastest rise time of 7 /rs 
was measured for TGF 101230.452 and gives an upper limit for the source radius of ~2km, 
assuming the speed of light [Briggs et al. , 2010] and 70 m assuming the propagation velocity of 
a lightning stepped leader. To explore the possibility that a TGF with a longer rise time could be 
modified by dead time or otherwise by the measurement process appear to have a 7 /rs rise time, 
we performed 1000 simulations of a TGF with cr^ set to 15//S (see § A2). Analysis of these 
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simulations show that the probability of observing a 7 //s rise time from such a TGF is <0.5%. 
Therefore, we conclude that it appears unlikely that TGF 101230.452 has an incident rise time 
which is substantially longer than the observed 7 /rs. 

5. Discussion 

Assuming that most TGFs are generated inside thunderclouds, there are two classes of TGF 
models that we consider: 1) lightning leader emissions, either in association with large scale 
thundercloud electric fields or fields generated by the lightning leader, and 2) the relativistic 
feedback discharge (RFD) model. The former involves energetic electron emission from the 
high fields associated with the lightning leader. Based upon this mechanism, Celestin and Pasko 
[2012] suggest that TGFs are the result of either single or multiple overlapping pulses, with 
durations of about 10 /as, that are smeared out by Compton scattering in the atmosphere. On the 
other hand, for the RFD model [Dwyer, 2012], the production of relativistic runaway electron 
avalanches becomes self-sustaining through the generation of backward propagating runaway 
positrons and back-scattered x-rays. 

The TGFs detected by GBM exhibit varied temporal behavior and 19% of the sample consist 
of multiple pulses, which may be overlapping or isolated by up to several ms in some cases. 
The overlapping pulses seen in 31 TGFs in the sample are evidence that the TGF production 
process is capable of producing multiple pulses on very short timescales. The shortest resolvable 
peak separation time observed for overlapping GBM pulses is ~ 200/as. Multiple TGF pulses 
may be signatures of multiple lightning discharges. Simulations of the RFD model show that 
multiple pulses are routinely produced but requires a large potential difference in the avalanche 
region [Dwyer, 2012]. Overlapping pulses caused by longer feedback times in a pulse train are 
also predicted by the model. Due to the larger potential required for multiple events, they are 
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predicted to be rare. However, this GBM sample is comprised of those TGFs which triggered 
the GBM detector over a timescale of 16 ms and so is biased towards multiple events. 

Using simulations we have shown that the temporal parameters are relatively unaffected by 
dead time and pulse pile-up (Appendix A2). The duration derived from the pulse parameters 
should be considered an upper limit, with the true duration being ~ 10% shorter. With these 
caveats in place, the temporal parameters presented here can be used to test theoretical models. 
For a model in which the applicable velocity is the advancement of a lighting leader, the source 
radius in constrained to be < 100 m (Table 1). As this limit is comparable to the size of a 
stepped leader, it is easily compatible with lightning leader models for TGFs. For the model in 
which TGFs are the result of lightning leaders propagating into a high field region and initiating 
RREA multiplication [Dwyer, 2008], the size limit suggests that only one lightning channel is 
involved in the gamma-ray emission. However, as TGFs propagate at the speed of light, this 
constraint is weakened if several leaders emit simultaneously (closer in time than the minimum 
observed variability time-scale). Whether or not multiple channels would be distinguishable 
depends on the relative geometry of the branches and the location of Fermi. For the fastest rise 
times observed by GBM, multiple channels would be indistinguishable if the paths from the 
sources to GBM differ in length by < 2km. 

Dwyer [2012] modelled the relativistic feedback processes self-consistently and found that, 
for the cases investigated, longer TGFs often make fewer runaway electrons than shorter pulses. 
The rise time at spacecraft altitude is less affected by Compton scattering than the duration and is 
relatively unaffected by dead time (Appendix A2). Therefore in Figure 6 the rise time is plotted 
as a function of the observed pulse fluence. There is a slight tendency for pulses with longer rise 
times to have lower fluences, in qualitative agreement with the RFD model. However, pulses 
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with shorter rise times have fluences which span the observed range. Drawing conclusions from 
this is difficult, as the observed fluence is modified both by dead time losses and the unknown 
GBM-source geometry. The range of observed durations of TGFs has been attributed to short 
overlapping emissions at the source, which are then smeared out by Compton scattering [Ce- 
lestin and Pasko, 2012]. Any such model would have to explain the symmetric pulses observed 
by GBM and the significantly shorter rise times observed for asymmetric pulses. 

The RFD model predicts that the time profile is always roughly symmetric at the source and at 
higher energies (> 1 MeV) at spacecraft altitude. At lower energies the propagation of photons 
through the atmosphere will cause an elongation in time due to Compton scattering. The rise 
times of the asymmetric pulses are significantly shorter than the symmetric pulses while their 
fall times are comparable. The observed rise and fall times of GBM TGFs can be explained 
by the RFD model if, at the source, pulses are intrinsically symmetric and have a range of 
durations. In this scenario short TGFs have tail emission which is easily observed and are 
detected as asymmetric pulses with short rise times. The Compton tails of longer pulses may be 
masked by the main pulse emission, particularly for weak TGFs with fewer statistics, and are 
observed as symmetric pulses with longer rise times. 

The fall times observed range from ~50 - 200yus. The observed fall time is the combination 
of the intrinsic fall time of the pulse and an additional Compton scattering component. The tail 
emission is also observed to be longer due to dead time since the peak emission is suppressed 
more than the tail, which increases the observed fall time of asymmetric TGFs by ~10% (Fig¬ 
ures Appendix A2). For a given range of altitudes and source - detector geometries, 0stgaard 
et al. [2008] showed that time delays on the order of 100//s are expected from Compton scat- 
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tering. With the inclusion of instrumental effects, the observed asymmetry can therefore be 
explained by Compton scattering. 
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Appendix A: Appendix 

We use Monte Carlo simulations to address issues in the data that require further examination. 
Firstly, the degree to which the observed parameters are affected by dead time must be quanti¬ 
fied. We also test if pulses with longer rise times than fall times can be explained by statistical 
fluctuations. These are addressed in turn below. 

Al. Simulation method 

For the purposes of simulation, the TGFs were divided into background and pulse compo¬ 
nents. The background was simulated as a homogeneous Poisson process, with a rate which 
was set to the mean of the observed background rate distribution for the Nal and BGO detec- 
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tors, which is 1.1 and 1.6 KHz respectively. Eaeh eount was assigned a eorresponding energy 
value whieh was randomly drawn from the observed baekground eount speetrum distribution. 

The pulse eomponent was simulated by first assuming a partieular pulse profile and then using 
the thinning method for a non-homogeneous Poisson proeess {Lewis et al, 1978]. Eaeh pulse 
eount was assigned a eorresponding energy from a speetrum whieh was derived by propagating 
the photons resulting from the RREA of an instantaneous distribution of eleetrons in an eleetrie 
field of 400 kV/m from an altitude of 15 km to 565 km in a 45° beam. This simulated photons 
were integrated in an annulus of 25 km at a sub-satellite to souree distanee of 200 km, eorre¬ 
sponding to the maximum probability of TGE deteetion (Briggs et al. [2013] and referenees 
therein). The Nal energy values were drawn in the range 4 keV - oo and the BGO energy values 
were drawn from 100 keV - oo. The pulse and baekground eomponents were then eombined. 
The eombined time profile eould then be passed through a filter whieh emulates the deteetion 
proeess. This filter ineludes effeets from both energy-dependant dead time and pulse pile-up 
effeets. 

A2. Effect of dead time on observed parameters 

Simulations were performed of symmetrie and asymmetrie time-profile TGEs with different 
durations and pulse shapes. Eor eaeh set of temporal properties, two different peak rates were 
simulated, one representative of the population and one typieal of the brightest events. These are 
referred to as ‘representative’ and ‘bright’ in the text. In total, 8 pseudo TGEs were simulated. 

Eor the symmetrie profiles, the parameters were set to crr=(Tf=l00ns and crr=(Tf=200/us 
for peak rates of 75 and 100 KHz. Eor the asymmetrie profiles, the parameters were set to 
cr,.=25 /IS, (Tf =125 /IS and (rr=50/is, af=l00/is for peak rates of 150 and 200 KHz. Eaeh TGE 
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was simulated and passed through the instrumental filter 1000 times. The resulting TGFs were 
pulse fit and the distribution of observed parameters compared to the input values. 

The observed peak rate and count fluences were heavily modified for both symmetric and 
asymmetric TGFs, with the observed value typically differing from the input value by ~ 10 - 
30 % depending on both the duration and initial rate. Due to their greater rates, the effect is 
more pronounced for the asymmetric TGFs. 

For the temporal parameters, the means of the measured distributions are broadly consistent 
with the input values, with a tendency to shift to longer times. This shift is generally small, with 
the input value usually contained within +1 cr of the mean. The exception is the fall times for 
the asymmetric TGFs, which show the same trend but to a greater degree. This is illustrated 
in Figures 7 and 8 which show the observed rise and fall time distributions respectively for a 
strong asymmetric TGF with the input value represented by a vertical black line. 

This elongation has the effect of artificially increasing the observed duration, however, the 
effect is small, even for very asymmetric high-rate TGFs, with the mean of the observed du¬ 
ration distribution typically differing from the input value by 20 - 30 fis. This is highlighted 
in Figure 9, which shows the observed distribution for an bright TGF with an asymmetric time 
profile. These results are consistent with that of Tierney et al. [2013], which showed that the 
observed temporal properties were not significantly effected by dead time and pulse pile-up. 

Another point of interest is the extent to which dead time can change a pulse from symmetric 
to asymmetric and vice-versa. For the asymmetric TGFs that were strongly asymmetric (Sxcr,. = 
(Tf) the effect was very small, with less than 1 % changing from asymmetric to symmetric. For 
those which were only slightly asymmetric (2 x cr,. = cry) the effect was more significant, 
with approximately 2 % of the representative TGFs and 5 % of the bright TGFs changing from 
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asymmetric to symmetric. For the symmetric TGFs, the discussion is slightly more nuanced. 
Due to statistical fluctuations (discussed in more detail in the following section), only 75% of 
the sample before dead time is applied are consistent with symmetry. Once dead time filter is 
applied and the TGFs refit, this increases to 80%. 

Several of the observed TGFs exhibit extremely fast temporal variation, with rise times less 
than 10 fis. For the simulations discussed above, the incident and observed rise times were found 
to be in reasonable agreement. However, it is possible is that a bias exists for shorter timescales 
and these short events are due to TGFs with longer rise times that have been modified by dead 
time. To explore this possibility, we consider the TGF with the shortest rise time, 101230.452, 
which has an observed CTr and cr/ of 6 and 71 //s respectively. Scaling the peak rate to 250 kHz 
and setting cr,. to 15/iS, we simulate and pulse fit this profile 1000 times. The distribution of 
rise times are consistent with a Gaussian. Based on the fit parameters, we can conclude that the 
probability of obtaining a 7 rise time from the distribution is ~0.5%. The distribution of rise 
times and the Gaussian fit are shown in Figure 10. 

A3. Significance of pulses with longer rises than falls 

Approximately 5 % of the pulse fits to the data result in longer rise times than fall times, 
that are not consistent with symmetry within the rise and fall time uncertainties. To determine 
whether this could be explained by statistical fluctuations from the symmetric group, 20 sets of 
pulse parameters were randomly selected from the sub-sample of fits consistent with symmetry. 
Each of these TGFs were then simulated 100 times and pulse fit. Figure 11 shows the result¬ 
ing rise time vs fall time. Approximately 12% of the sample have longer rise time than fall 
time. The equivalent value for the data is ~ 15 %, mostly consistent with the simulations. This 
implies that the majority of the TGFs pulse fits with longer rise than fall times are consistent 
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with statistical fluctuations from the symmetric group. The remainder are likely unrecognised 

overlapping pulses that were fit with a single Norris pulse. 

Notes 

1. http://F^rm/.gsfc.nasa.gov/ssc/ 

2. http://gammaray.nsstc.nasa.gov/gbm/science/tgf 
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Figure 1. T 90 duration calculated from the count data versus that derived from the pulse fits. The 

values from the two methods are eonsistent. 


Table 1. Constraints on the souree emission radii for the 3 TGF pulses in the sample with the shortest 
rise times assuming different propagation velocities. 


TGF 

Rise Time (/rs) 

Propagation 


Maximum Souree Radius (m) 

TGF 101230.452 

7 + 3 

Speed of light (3 x 10^ m/s) 
Stepped lightning leaders (~ 

1 X 10’ m/s) 

2100 + 900 

70 + 30 

TGF 100716.443 

8 + 4 

Speed of light (3 x 10^ m/s) 
Stepped lightning leaders (~ 

1 X 10’ m/s) 

2400+ 1200 

80 + 40 

TGF 100426.657 

9 + 3 

Speed of light (3 x 10^ m/s) 
Stepped lightning leaders (~ 

1 X 10’ m/s) 

2700 + 900 

90 + 30 
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Figure 2. GBM TGFl 10828.435 consisting of multiple and overlapping pulses. The individual Norris 
pulse fits are shown as dotted lines; the summed fit is shown as a solid line. Sinee the fits are performed 
at the native resolution of 2 /us, the ehoiee of data binning for the lighteurve is for illustration purposes 
only and does not have any effeet on the pulse fit. 
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Figure 3. Rise time versus fall time for 250 pulses fit with Norris functions. The dashed line represents 
the line of equal rise times and fall times. 
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Figure 4. Rise time versus fall time for pulses fit with Gaussian funetions (purple) and pulses best fit 
with lognormal funetions (gray). The dashed line represents the line of equal rise times and fall times. 
The pulse parameters obtained from the simpler Gaussian and lognormal fits are better eonstrained than 
those derived from the Norris funetion. 
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250 TGF pulses. 
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Figure 6. Fluence versus rise time. The observed fluenee is a lower limit on the aetual value due to 
instrumental effeets. There is a slight tendeney for pulses with longer rise times to have lower fluenees, 
eonsistent with the RFD model. However, pulses with shorter rise times have fluenees whieh span the 
observed range. 
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Figure 7. Distribution of observed rise times for a strong (input peak rate of 200 KHz per detector) 
asymmetric TGF. The input simulation rise time is plotted as vertical black line. The observed times 
are completely consistent with the input. 
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Figure 8. Distribution of observed fall times for a strong (input peak rate of 200 KHz per deteetor) 
asymmetrie TGF. The input simulation fall time is plotted as vertieal blaek line. The observed times 
are signifieantly shifted to longer times. However, the effeet is relatively small at ~10 %. 
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Figure 9. Distribution of observed durations for a strong (input peak rate of 200 KHz per deteetor) 


asymmetrie TGF. The measured duration is clearly shifted to longer times, but the effect is small, with 


the mean being shifted by ~ 20 yU5. 
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Figure 10. Distribution of rise times for a strong TGF with a short rise time (cr^ =15 yus). The solid 
line is the result of a Gaussian fit to the data and the dashed vertieal line indieates the 7 /rs value. The 
probability of the 7 yus value being drawn from the distribution is ~0.5%. 
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Figure 11. Rise time vs fall time for 20 symmetrie TGFs, eaeh simulated 100 times. Points eonsis- 
tent with symmetry are shown in blue, those with longer rise times than fall times are shown in red. 
Approximately 12 % of the sample have longer rise time than fall time. 
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